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Abstract 



We investigated the influence of dark matter on light propagation in the solar sys- 
tem. We assumed the spherical symmetry of spacetime and derived the approximate 
solution of the Einstein equation, which consists of the gravitational attractions 
caused by the central celestial body, i.e. the Sun, and the dark matter surrounding 
it. We expressed the dark matter density in the solar system in the following simple 
power-law form, g{t,r) = p{t){l/r)^ , where t is the coordinate time; r, the radius 
from the central body; £, the normalizing factor; k, the exponent characterizing 
r-dependence of dark matter density; and p{t), the arbitrary function of time t. On 
the basis of the derived approximate solution, we focused on light propagation and 
obtained the additional corrections of the gravitational time delay and the relative 
frequency shift caused by the dark matter. As an application of our results, we 
considered the secular increase in the astronomical unit reported by Krasinsky and 
Brumberg (2004) and found that it was difficult to provide an explanation for the 
observed dA\J/dt = 15 ± 4 [m/century]. 
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1 Introduction 



The existenc e of dark matter was fi rst indicated bv IZwickvl (119331 ) and 
subsequently by Rubin and Ford ( 19701 ): Rubin et al. ( 1980 ). According to the 
recent cos mological observatio n, i.e. Wilkinson Microwave Anisotropy Probe 
(WMAP) (ISpergel et al.l . l2003l ) . it is suggested that the majority of mass in our 
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Universe is the dark matter, which is approximately 6 times that of ordinary 
(baryonic) matter. Furthermore, it is considered that dark matter has played 
an important role in the large-scale structure formation of the Universe. Actual 
observations such as 2dF[irand SDS^Zl and results of numerical simulations, 
which assume the existence of dark matter such as Virgo Consortiunf^. are in 
good agreement with each other. Therefore, dark matter is considered to be 
the fundamental material of our Universe, even though its details are presently 
cloaked in mystery. 



If dark matter is an essential component in the Universe, it is interesting and 
worthy to investigate its existence and detectability in our neighborhood area, 
the solar system. The number of recent astronomical and astrophysical mea- 
surements in the solar system have drastically increased, and they have led to 
(i) deep understanding of planetary dynamics and fundamental gravitational 
physics, (ii) significant improvement in lunar and planetary ephemerides, and 
(iii) precise determination of various astronomical constants. 



To date, the gravitational influence of dark matter on planetary motion, 
such as the additional perihelion advance, has been studied by several au- 
thors. For instance, the e f fect of galactic dark mat t er has been considered 



bvlBraginskv et all (19921). iKlioner and Soffell (11993f ). iNordtvedtl (119941 ). and 
Nordtvedt and Soffell (119951 ). On the other hand, the upper limit o f the dark 



matter density in the so 



ar system has been estimated within a range (jAnderson et al. 



1989; 


Gr0n and Solene. 


1996; 


Khriolovich and Pitieva 


2006; 


Khriplovich . 


2007 


; Frere et al.. 


2008 


) such that 



pir^ < lO-^*^ ~ 10-2° [g/cm% (1) 
Recent l y, the upper limit of the planet-bound dark matter was also evaluated 



(lAdleii . l2008a 



However, its contribution to light propagation has hardly been examined, in 
spite of the fact that the current accurate observations have been archived by 
an improvement in the observation of light/signal; round-trip time (radar/laser 
gauging techniques and spacecraft ranging), amelioration of atomic clocks, ra- 
dio links of spacecraft, and increasing stability of frequency standar d. More- 

over, planned space missions, such as GAIAH, SIM0, LISA0, LATOR (iTurvshev et al.l 



20041 ). and ASTROD/ASTROD-1 (InI. l2007h . require accurate light propaga- 
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tion models. Therefore, it is noteworthy to examine how dark matter in the 
solar system affects light propagation and whether its traces can be detected 
if it really exists in our solar system. 

On the other hand, because of a drastic improvement in the measurement 
techniques used in the solar system, som e unexplained p r oblern atic phenom- 



ena occur red, such as pio n eer a nomaly (j Anderson et al.l . Il998l ). Earth flyby 



anomaly (lAnderson et al.l. 120081) , secular increase in the astronomical unit 



(IKrasinsky and Brumbergl . |2004J ). and anomalous perihelion precession of Sat- 
urn flloriol . l2009l l Currently, the origins of these phenomena are far from clear, 
nevertheless, dark matter may cause some significant coii t ribut ion to these 



phenomena fiNietol . I2OO8I : lAnderson and Nietol . lioool : lAdlerl . l2009h 



In this study, we will examine the influence of dark matter on light propaga- 
tion in the solar system. First, we assume the spherical symmetry of spacetime 
and derive a simple approximate solution of the Einstein equation, which con- 
sists of gravitational attractions caused by the central celestial body, i.e, the 
Sun, and the dark matter surrounding it. Then, we will focus on formulating a 
light propagation model and estimating the additional effects of gravitational 
time delay and the relative frequency shift of a signal. As an application of 
our results, we will consi der the secular increase iri the as tronomical unit (of 
length), AU, reported by IKrasinsky and Brumbergl (120041 ). 



This paper is organized as follows. In Section [21 we explain the model of 
spacetime and some assumptions. In Section [31 we derive the approximate 
solution of the Einstein equation. In sections [H and \5[ we investigate time 
delay and relative frequency shift, respectively. In section [6l we focus on the 
application of our results to the secular increase in the astronomical unit. 
Finally, in section [71 we provide the summary of our study. 



2 Model and Assumptions 



Before deriving the approximate solution of the Einstein equation, we ex- 
plain the model of spacetime and some assumptions. First, we suppose that 
the spacetime is characterized by the gravitational attractions caused by the 
central celestial body, i.e. the Sun, and dark matter surrounding it. Then, we 
express the spherically symmetric form of metric as 

ds^ = -e^c^df + e^dr^ + r'^^dO'^ + sin^ ^d^^), (2) 

where and are functions of time t and radius r, and c is the speed of 
hght in vacuum. 
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As the stress-energy tensor TS" we presume the following form 



T° = -^(^,r)c^ T^ = a{t,r)c, T] = 0. (3) 

Tq is related to the dark matter density g{t, r). At this time, since we do not 
have any knowledge about the actual distribution of dark matter in the solar 
system, then we suppose the following simple power-law form: 




Q{t,r)=p{t)\-\ , (4) 



where ^ is the normalizing factor that is chosen as £ = and is the orbital 
radius of Earth, k is the exponent characterizing r-dependence of Q{t, r), and 
p{t) represents the time variation of the dark matter densit}0. The test par- 
ticle, i.e. photon, is subjected to the gravitational attractions caused by the 
central body and dark matter, which is confined within a spherical shell of 
radius r at time t. p{t) is an arbitrary function of time t, however, for the sake 
of simplicity, we assume that the time variation of the dark matter density in 
the solar system is considerably slow; 



- PC + I 



(«-«o), (5) 



where subscript denotes the initial epoch of planetary ephemerides. As the 
one possibility, the time variation of the dark matter density g{t, r) may be 
caused by the motion of the solar system in our Galaxy, if the distribution of 
galactic dark matter is inhomogeneous. 

The dark matter density g{t, r) observed in the solar system may be gener- 
ally expressed as a sum of solar system-bound (or local) dark matter §^^"^^"^^1, r) 
and galactic dark matter r) as follows: 

Q{t,r) = 0^'°^'''^ {t,r) + ^(g^i^^"^) (t, r) . (6) 



Here, to simplify the situation, we assume that the spacetime is spherically 
symmetric and the time variation of the dark matter density g{t, r) is caused by 
the inhomogeneity of the galactic dark matter as mentioned above. Therefore 
it is possible to express 



' In this paper, Greek indexes run from to 3, and Latin ones do from 1 to 3. 
^ pit) can be considered to be the dark matter density observed around the Earth's 
orbit since £ = rE- 
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Q 



(galactic) ^-j 



(galactic) , dp 

Po + — 




In this case, po and dp/dt\o in ([5]) are 



(t - to) 



po = p(^°i^'^)+pfe^i^^"^), 



(ip 
'di 



dp 



(galactic) 



dt 



(7) 



(9) 



According to the recent investigation, i.e. iBertone and MerrittI (120051 ). the 
galactic dark matter density is of the order of 10~^^ [g/cm^]. This is several 
orders of magnitude smaller than the evaluated density of dark matter in the 
solar system (see ([T])). Therefore in this study, we suppose that 



(solar) dp 
P' = P^ ■ -dt 



dp 



(galactic) 



dt 



(10) 



Tq,T^ represent the time variation of energy and momentum flux. We auxiliary 
introduced these components to preserve the time dependency of the obtained 
solution. Though represents the stress part attributed to the dark matter, 
currently, the equation of state of dark matter is not known, therefore, we 
adopt the standard assumption that the dark matter is pressure-less dust 
particles p ~ and that its time variation is also negligible dp/dt ~ 0. 

Finally, we consider the choice of exponent k of g{t,r). Because the distri- 
bution of dark matter in the solar system is poorly understood, we adopt the 
following three indexes as examples: k = 1 (density decreasing with r), k = 
(constant density), and k = —1 (density increasing with r). The density g{t, r) 
should be damped at a certain radius from the Sun (especially when k = 
and k = —1), and it is natural to imagine that g{t,r) reaches asymptotically 
for the galactic dark matter density, i.e. ~ 10~^^ [g/cm^]. However, we are now 
interested in the astronomical observations within the quite inner (planetary) 
area of the solar system (see Fig. [T] for the conceptual diagram). Therefore, in 
this study, we do not consider the details of dark matter density distribution 
far away from the Sun. 



3 Approximate Solution of the Einstein Equation 



On the basis of the assumptions in Section [21 we obtain the approximate 
solution of the Einstein equation 
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Fig. 1. Conceptual diagram of r-dependence of dark matter density for A; = 1, A; = 0, 
and k = —\. We simply assume that the actual observations are carried out within 
the quite inner (planetary) area of the solar system <C r^. 



SvrG 



:ii) 



where G is the Newtonian gravitational constant, is the Einstein tensor, i?^ 
is the Ricci tensor, R is the Ricci scalar, and is the stress-energy tensor in 
([3]). The non-zero components of the Einstein tensor are expressed as follows: 



G\- 



+ 



Gl 



e~ 



ldv\ 
r dr J 



n 



1 



2c2 



d^u 1 



\ dr ^ 
dt 



1 du dfi ^ 1 

2 dr dr r 



1 (9/i dv 
2~dt~dt 



cr dt^ ° cr dt 



'(9/i 
dr 



du^ 
dr , 



[12) 



(13) 



(14) 
(15) 



From the 00 component of the Einstein equation, we have 



r 1 - e-^) = m{t) + , , r^'", (16) 

where m{t) is an arbitrary function of time t, however, we choose m{t) such 
that it reduces to the Schwarzschild radius m{t) — m = 2GM/c^ = constant 
when g{t, r) = 0. Therefore, we obtain 
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1 - 



2GM MGp{t)t 2_fc 



3- A; 



(17) 



Using the 00 and 11 components of the Einstein equation, it follows 



(1^ 



where we kept the 0{c^'^) order terms only on the right-hand side. We obtain 
the following equation by integrating (1181) with respect to r, combing it with 
(1171) . and omitting the (9(c~^) and higher order terms: 



fit) 



2GM ^ 8nG p{t)t 



„2-fc 



(19) 



c2 {2-k){?>-k) 

Although fit) is also an arbitrary function of time t, we replace the time 



coordinate with J f{t)dt dt and delete f(t). Finally, we obtain 



2GM SnG 

+ 



2-k \ 2 7.2 



1 - 



c2 {2-k){?,-k) 
2GM 87TGp{t)i 



c'dt' 



c^r 



c2 3 -A; 



where dVt^ = dO"^ + sin^ 6*^02. When = and p it) 



Po 



(20) 



constant, (!20!) 



reduces to the metric derived by lGr0n and Solengl (119961 ): therefore, our solu- 
tion (l20l) is considered to be some extension of the solution obtained by Gr0n 
and Soleng. In the case of a solar system experiment, it is sufficient to use the 
following approximate form: 



ds^- 
where 



l--f/(t,r) 



l + ^Vit,r] 



dr^ + r^dVL^ 



(21) 



U{t,T) 
Vit.r) 



GM A'KGpit)^ 



,2-fe 



r (2-fc)(3-A;) 
GM A-KGpit)t 2_fc 



+ 



3 - A; 



(22) 
(23) 



We mention here that it is easy to incorporate the cosmological constant A 
in (l20l) or (12T!) . However, to focus on the effect of dark matter, we omit the 
— Ar2/3 term. 
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Fig. 2. Light/signal path. We assume that the first approximation of the light path 
is rectilinear along the x-direction (bold dashed line), b is the impact parameter, 
and r = \J + 6^. The actual light path is drawn by the bold solid line. 

4 Gravitational Time Delay 

J^.l Time Delay in Coordinate Time 



In this section, we calculate the time delay attributed to the dark matter. In 
the static spacetime, we can easily relate an affine parameter A to coordinate 
time t using by the Euler-Lagrange equation of (?00) i-e. Chapter 8 of [Weinberg 
( I1972I ). However, (l20l) or (12T!) is the non-static or time-dependent then it is 
not easy to calculate the geodesic equation analytically in general. Therefore, 
we consider an alternative approach. 

To begin with, we transform (1211) from spherical coordinates to rectangular 
coordinates. By the usual coordinate transformation, i.e. 



X = r sin Q sin 0, y = t sin 9 cos 1 



r cos 6, 



(24) 



(1211) is rewritten as (e.g.. iBrumberg] (jl99ll )) 



(25) 



where 5ij is the Kronecker's delta symbol. We suppose that the actual light 
path is calculated along with the approximate rectilinear path (x-direction) 
such that 



y = b = constant, z = 0, r = 



(26) 



where b is an impact factor (see Fig. [2]). Hence, (|25|) becomes 



ds^ = -il- lu] c^de + f 1 + ^ dx\ 



(27) 
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The world line of the light ray is null geodesic ds^ = 0, therefore from (127|) we 
have 



dt 
dx 



1 + 



GM 



1 + 



{2-k){3-k) 



3-k 



i.28) 



We express p{t) in the form of 



To obtain the round-trip time from (l28i) . we assume that the lapse time At 
is expressed by a linear combination of each effect as follows: 



At 



AtpN + At + At^m 



(29) 



where AtpN corresponds to the Shapiro time delay in 1st post-Newtonian 



approximation, At 



(const) 
dm 



is attributed to the static part of the dark matter 



density (po of ([5])), and At^^l^ is the contribution of the time-dependent part 
of dark matter density {dp/dt\ot of ([5])). The post-Newtonian parts are easily 
obtained as follows: 



At 



X2 - Xi 



pN 



+ 



C 

GM 



2 In 




Xi 



xj + b'^ ^ 



(30) 



Next, we calculate the time delay caused by dark matter. The static part 
Atj™'^''*'* is straightforwardly integrated as 



At 



(const) 
dm 



^^^^H{xi,X2; k) 



.2ib' In f£±v^ 
\{xl — xf) — b'^{x2 — Xi] 



H{xi,X2] k) = { 



36^ In 



X2 + ^/xl+^P 
■3 



2(x2Vxi + 62 -x^Jxl + b'^ 



+ 3b'^{x2Jxl + 62 - xiJxl + 62) 



ik = l) 
{k = 0) 



{k = -l). 



(31) 



(32) 



Finally, we compute the time-dependent part At^j^j {dp/dt\ot part). If light is 
emitted from Earth at t = T and it reaches the reflector (planet /spacecraft) 
at t = T -|- At[i*j|j, then At^^l^ satisfies 
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T+At' 



(t) 



-dt = cln 
t 



dm _ '^P 

dt 



^2 

y [■ ■ ■] t^a^, 



(33) 



where the integral [■ ■ ■] dx is the same as that in the case of the static part 
Atfr'*^ Since dp/dt\^ < 1, Atf^ can be expressed as 



At 



{t) _l dp 



X2 

T j [■■■] dx. 



<^"^ c dt 

Therefore, the time delay caused by dark matter is expressed as 



(34) 



Atfjm = At 



(const) 
dm 



dp 



T H{xi, X2] k). 



(35) 



In this case, the time delay caused by dark matter can be characterized by 
the density at the emission time of signal t = T, that is p(T) = po + dp/ dt\oT . 

In order to apply (|5U]) and (jHS]) to two-way light propagation, i.e. Earth 
— i> receiver (planet /spacecraft) — Earth, we consider the following situation; 
the light path used in computation is parallel to the x-axis. Earth, and the 
receiver (planet /spacecraft) are located x = oe and x = — or, respectively 
(see Fig. [2] again). We suppose that (a) during the round-trip of light. Earth 
and the receiver are almost at rest and that (b) the time variation of the dark 
matter density is also considerably slow. In other words, the time lapse AT 
can be mainly determined using by the dark matter density at the emission 
time t = T, p{T). Hence, the round-trip time in the coordinate time, AT is 
expressed as 



, ^ aE + Or 2GM 
AT = 2— H 



2 In 



«E + V«E + b^) (or + V + 

62 




27rG f dp 



T H{aE, Or; k), 



(36) 



where we substitute 7i(aE,ciR;fc) = H{0,aE',k) + iJ(0, or; /c). To calculate 
([36|). we refe r red to an approach in Section 40.4 and Figure 40.3 shown by 
Misner et aP fll970h . 



Let us estimate the order of the time delay ATdm- Because dp/dtloT is now 
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1e-14 


yfc = -1 




1e-16 




k = 


1e-18 






, CO , 






1 1e-20 

< 

1e-22 






1e-24 






1e-26 





5 10 15 20 25 30 35 40 45 50 

an [AU] 

Fig. 3. aR dependence of additional time delay ATdm- As po, we adopt 
po ~ 10-^*^ [g/cm^] and fixed ue = 1.0 [AU](= 1.5 x 10^^ [m]) and 
b = 0.001 [AU](= 1.5 X 10^ [m]). 

anticipated to be considerably smaller than the dominant part po, we neglect 
the dp/dt\oT term here and evaluate 

ATarn^^^n{aE,an;k). (37) 

Fig. [3] illustrates the or dependence of the time delay ATdm- We adopted po ~ 
10~i^ [g/cm^], which is the largest upper limit obtained from the dynamical 
perturbation on planetary motion. We fixed oe = 1-0 [AU](= 1.5 x 10^^ [m]) 
(orbital radius of the Earth) and impact parameter b = 0.001 [AU](= 1.5 x 
10^ [in]). If the dark matter is accumulated in the neighborhood of the Sun 
{k = 1), ATdm ~ 10~^^ [s] in a given range of or. When k = and k = — 1, 
ATdm is of the order of 10~^° [s] in the inner planetary region, while in the 
outer planetary region, it is of the order of 10"^^ < ATdm < 10^^'' [s] (A; = 0) 
and 10^^^ < ATdm < 10^^^ [s] {k = —1). However, the current observational 
limit in the solar system is ~ 10~^[s] or a few 100 [m] for planetary radar and 
10~^^[s] or a few [m] for spacecraft ranging; the internal error of the atomic 
clocks on Earth is ~ 10~^ [s]. Then, at this time, it is difficult to extract the 
trace of dark matter from the ranging data. 

4-2 Time Delay in Proper Time 

The round-trip time of the light ray (136|) is expressed in the coordinate 
time. However, the actual measurement is performed by the atomic clocks on 
the surface of Earth, which shows proper time r. Therefore, we must transform 
(I3HD into proper time. Presently, it is sufficient to use the equation of proper 
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time for the quasi-Newtonian approximation such that 



dr 



1 



— = l--{U+-v 
at c 



(3J 



Evaluating dr/dt around the orbit of Earth and keeping the 0{c ^) terms 
only, the round-trip time Ar measured in proper time is given by 



dr 

At= — 

dt 



AT 



oe + ciR 2GM 

2 ~t~ Q 



Oe 



21n- 



Or 



flE + V4 + &^)(«R + V«R + 



62 



-2 



a| + 62 ^ V4, + &^ 

oe + Or 1 9 GM 
2 oe 



+ 



27rG' 



dp 



c3 yP'^Tt 



T 7^(aE,aR; A;) 
/ 



T 



(39) 



where we is the orbital velocity of Earth. 



5 Relative Frequency Shift 



We use 
defined as 



to derive the relative frequency shift of signal ?j_J, which is 



5u 

V 



dAT 
dt 



(40) 



When the light ray pa s ses ne ar the limb of the Sun such as in the Cassini exper- 
iment teertotti et allboO^h . the conditions qe, or ^ b, da-^/dt, andda-Ei/dt <^ 



db/dt hold, where b = yfeo + ("^i^)^- Then, the relative frequency shift caused 
by the Sun, ?/pN, and dark matter, ?/dm, are expressed as 



y = ypN + Vdn 

_ 8GM db 



(41) 
(42) 



^ Here, y is not th e y-coordinate, but the relative frequency shift according to 
Bertotti et all (|2003l ). 



12 



1e-25 


-1e-25 
-2e-25 



k = 
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-5e-25 
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an [AU] 

Fig. 4. Additional relative frequency shift caused by dark matter, and ydm is plotted 
as a function of or. We set bo 
t = l [day]. 



/C(aE,aR; /c) 



8bi (in 



4aEQR 
b2 



2i?sun,i?Sun ^ 6.9 X 10« H>^^ - 30 [km/s], and 

(43) 



1 ) - 



|6(aE + aR)| 



ik = l) 
{k = 0) 



(44) 



1^ 



6^ In ^ - (a| + a|)l f (k 



Fig. m shows the relative frequency shift caused by dark matter, j/dm as a func- 
tion of Or. In this plot, we substitute 60 = 2i?sun5-Rsun — 6.9 x 10® [m],t> ~ 
30 [km/s], and t = 1 [day] [t = gives the closest point). The order of magni- 
tude of ?/dm is ~ 10~^^; however, currently the stability of frequency standard 
is of the order of 10~^^ or even higher. Therefore, the expected frequency shift 
caused by dark matter is approximately 10 orders of magnitude smaller than 
the present observational limit of frequency. 



6 Application to Secular Increase in Astronomical Unit 



In this section, we apply the previous re sults to the secular increas e in th e 
astronomical unit (of length) reported by iKrasinsky and BrumbergI (120041 ). 
The astronomical unit (AU) is one of the important scales in astronomy, and 
it is the basis of the cosmological distance ladder. AU is also a fundamental 
astronomical constant, which gives the relation between two length units; 1 
[AU] in the astronomical system of units and 1 [m] in SI ones. Presently, AU is 
determined by using the planetary radar and spacecraft ranging data (round- 
trip time of light ray), and the latest best-fit value is obtained as (jPitjeval . 
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20051 ) 



1 [AU] = 1.495978706960 x lO" ± 0.1 [m]. (45) 

We use the calculated planetary ephemerides (solution of equation of motion) 
to compute the theoretical value of the round-trip time ttheo using the following 
formula: 



ttheo = ^AU[s], (46) 
c 

where citheo [AU] is the interplanetary distance, ttheo [s] is compared with the 
observed round-trip time tobs [s], and AU is optimized by the least square 
method. 

However, when Krasinsky and Brumberg replaced ttheo with 



^theo 



''theo 



AU + ^(t-to) 



(47) 



and fitted it to the observational data, they found that dAU/dt had a non-zero 
and positive secular value, 15±4[m/century], where to is the initial epoch. The 
evaluated value dAU/dt = 15±4 [m/century] is approximately 100 times that 
of the current determination error of AU (see (l45l) ). At present, the time de- 
pendent part {dAU /dt){t — to) cannot be related to any theoretical predictions, 
hence, several attempts have been made to explain this secular increase in AU 
on the basis of various factors such as the effects of cosmological exp ansion 
( Krasinskv and Brumberg . 2004 : Mashhoon et al. . 2007 : Arakidal 2009 ). mass 
loss of the Sun (IKrasinsky and Brumberg . l2004l: iNoerdlingen. i2008[). a n d tirn e 



variation of the gravitational constant G (IKrasinsky and Brumbergl . |2004| ). 
However, unfortunately, thus far, none of these factors seem to be responsible 
for the secular increase in AU. 



It is noteworthy that the observed dAU/dt does not imply the expansion 
of planetary orbit and/or an increase in the orbital period of a planet. As a 
matter of fact, the determination error of the latest planetary ephemerides 
isconsiderably smaller than the reported dAU/dt (see in Table 4. of iPitjeva 
(120051 )). Hence, dAU/dt may be caused by some effects on light propagation, 
and not by the dynamical perturbation on planetary motion. 

Moreover, AU denotes not only the conversion constant of the length unit 
but also the value that characterizes the GM of the Sun in SI units such that 



GMsnn = k'AU'/d' [mVs'], 



(4J 
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where k = 0.01720209895 is the Gaussian gravitational constant, and d is a 
day such that d = 86400 [s]. Therefore, the observed dAU/dt may be related to 
an increase in the dark matter density such that GM(t) = G(Msun + ^dm(^)), 
where Mdm(t) is the total mass of dark matter within a planetary orbit at 
time t. 

Then, let us evaluate the extent of time variation of dark matter density 
dp/dt\Q in fl36|) that is needed to explain the observed dAU/dt. We have 



dtheo dAU 
c dt 



T 



27iG dp 
~lf 'dt 



TR^, 



(49) 



where we set Ti.{a-E, Qr; k) ~ R'^, and R is the orbital radius of a planet. In the 
case of Earth-Mars ranging, we let i? ~ 1.52 [AU] (orbital radius of Mars). To 
obtain the reported dAU/dt, dp/dt\o must be of the order of 10~^ [g/(cm^s)] 
and dp/dt\oT ~ l[g/cm^] for T ~ 100 [y]. However, this value corresponds 
to the density of water; therefore, this possibility of achieving such value is 
unrealistic and should be made an exception. 



7 Summary 



We investigated the influence of dark matter on light propagation in the 
solar system. We used the simplified model to derive the approximate solution 
of the Einstein equation, which consists of the gravitational attractions caused 
by the central celestial body, i.e. the Sun, and dark matter surrounding it. We 
found that the de rived metric (l?ID can be considered to be an extension of the 



previous work by lGr0n and SolengI (119961 ). We assumed that the simple time 



variation of dark matter density, and focused our discussion on light propa- 
gation then computed the additional corrections of gravitational time delay 
and relative frequency shift. However, the expected effects were considerably 
smaller than the current observational limits, even when we considered the 
largest upper limit evaluated from the planetary perturbation caused by dark 
matter, po ~ 10^"*^^ [g/cm'^]. 

We applied th e obtained results to the se c ular i ncrease in the astronomical 



unit reported by iKrasinsky and BrumbergI (120041 ) and considered the possi- 



bility of explaining the observed dAU/dt = 15 ± 4 [m/century] on the basis 
of the time variation of the dark matter density. We found that to induce the 
obtained dAU/dt, the change in the dark matter density dp/dt\o in ([5]) must be 
of the order of 10~^ [g/(cm^s)] and that dp/dt\oT ~ 1 [g/cm^] for the interval 
T ~ 100 [y]. However, it is completely unrealistic to achieve these values, and 
the existence of dark matter and its time variation cannot explain dAU/dt. 
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As mentioned in the previous section, some attempts were made to show 
the secular increase in AU. However, the origin of dAU/dt is presently far from 
clear. As the one possibility, it is believed that the most plausible reason for the 
origin of dAU/dt is the lack of calibrations of internal delays of radio signals 
within spacecrafts. Nevertheless, as other unexplained anomalies discovered in 
the solar system, dAU/dt may be attributed to the fundamental property of 
gravity, therefore, this issue should be explored in terms of all possibilities. 

Though it is currently impossible to detect the evidence of dark matter 
from light propagation, some planned space missions, especially ASTROD, 
are aimed to achieve a clock stability of 10"^'' over a travel time of 1000 [s] 



( iNil . 120071 ). Improvement in both the laser ranging technique and the clock 
stability may enable us to observe the trace of dark matter, if it really exists 
in the solar system. For this purpose, it is very important subjects to develop 
a rigorous hght propagation model. 

In particular, since it is not easy to analytically calculate the time-dependent 
null geodesic equation, in this study, we integrated fl28|) assuming the simple 
linear combinations of each effect. However, from the theoretical point of view 
and some astronomical and astrophysical applications such as formulation 
of the cosmological gravitational lensing in the expanding background, it is 
noteworthy to develop a method to analytically compute the time-dependent 
geodesic equation. 
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